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^ ; Abstract 

' ^ , We prove that adiabatic computation is equivalent to standard quantum computation even 

| ' when the adiabatic quantum system is restricted to be a set of particles on a one-dimensional 

chain. We give a construction that uses a 2-local Hamiltonian on nearest neighbors using 
particles that can have ten distinct states. This implies a construction of a one-dimensional 
chain of qubits in which the Hamiltonian is 6-local. We adapt this construction to show that 
the 2-local Hamiltonian for 13-state particles is QMA-complete which in turn implies that the 
8-local Hamiltonian restricted to a one-dimensional chain of qubits is QMA-complete. 
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53 ■ 1 Introduction 
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Adiabatic computation was introduced by Farhi, Goldstone, Gutmann and Sipser [?] as a means 
of solving difficult classical optimization problems. Although it appears likely that algorithms for 
this model require exponential time to solve NP-hard problems, adiabatic computation remains an 
■ appealing alternative to standard computation since it may be more robust against certain types 

^sO . of quantum errors [?]. 

The idea behind adiabatic computation is to define two Hamiltonians Hi n u and Hfi na [. The 
ground state of B.i n a (the eigenstate coresponding to the smallest eigenvalue) should be easy to 
prepare, like a tensor product states. The desired output is the ground state of Hj ina i which should 
'. somehow encode the solution to the problem. We also require that both Hamiltonians be local in 

that they are the sums of Hamiltonians, each of which operates on a constant number of particles. 
Thus, they are efficiently described by enumerating the matrices associated with each component. 
The Hamiltonian is slowly varied from Hi n a to Hfi na i by increasing s from to 1 in the following 
expression: H(s) = (1 — s)Hi n a + sHfi na i. The adiabatic theorem (expressed more formally in the 
next section) states that if a system begins in the ground state of Hj n # and the system is varied 
slowly enough, it will remain in the ground state. Thus, the final state will be the desired ground 
state of Hfi na i. The required time to vary s from to 1 depends on the minimum spectral gap of 
H(s). In order to think of adiabatic quantum computation as computing a classical function (as 
opposed to quantum states), a measurement of one or more of the particles is then performed to 
yield a classical output. 

van Dam, Mosca and Vazirani proved that the standard model of quantum computation is at 
least as strong as the adiabatic quantum model [?]. Subsequently, Aharonov et al. showed that 
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an arbitrary quantum circuit can be simulated by an adiabatic computation [?], thus establishing 
polynomial equivalence of the two models. Furthermore, the construction in [?] showed that any 
quantum computation could be efficiently simulated by an adiabatic computation with 2-local 
nearest-neighbor Hamiltonians operating on six-state particles on a two dimensional grid. This was 
later improved by Oliveira and Tehral [?] to a two-dimensional grid of qubits (two state particles). 
We prove the following theorem for 1-dimensional systems: 

Theorem 1 An arbitrary quantum circuit can be simulated by an adiabatic computation with 2- 
local nearest neighbor Hamiltonians operating on a 1-dimensional chain of ten-state particles. This 
in turn implies that the construction will work on a 1-dimensional chain of quibits with 6-local 
Hamiltonians, each of which operates on continguious sets of qubits. 

A slight modification of our construction has implications for the study of the complexity of 
quantum systems. The study of quantum systems from the point of view of their computational 
complexity is an important part of the more general research program aimed at understanding 
the power of quantum computing. In particular, many one-dimensional quantum systems appear 
to be more tractable in some ways than their two-dimensional counterparts [?]. This paper seeks 
to answer whether certain hardness results for two-dimensional quantum systems can be carried 
over to the one-dimensional case. In particular, can a one-dimensional quantum system have a 
sufficiently rich structure to implement universal computation? 

Proving that a given question about a quantum system is QMA-complete provides strong ev- 
idence that it is computationally difficult, even with the benefit of a quantum computer. The 
class QMA is the quantum analog of NP and MA. That is, it is the set of all languages that can 
be probabilistically verified by a quantum verifier in polynomial time. Kitaev initiated the study 
of QMA-complete problems by defining the local hamiltonian problem, the quantum analog of 
SAT [?]. In this problem, one is given a Hamiltonian and a guarantee that the lowest eigenvalue 
for H is either greater than some value b or less than another values a, where b — a is at least 
an inverse polynomial in the number of qubits in the system. The output of the problem is to 
determine which alternative is the case for H. Kitaev gave the first QMA-complete problem by 
showing that the 5-local Hamiltonian is QMA-complete [?]. This was later improved by Kempe 
and Regev who showed that the 3-local Hamiltonian is QMA-complete [?] and then by Kempe, 
Kitaev and Regev who showed that the 2-local Hamiltonian is QMA-complete [?]. Finally Oliveira 
and Tehral [?] showed that this latter result holds even when the 2-local Hamiltonian is constrained 
to be nearest-neighor interactions on a grid of qubits. They conjecture that it is not possible to 
extend this result for a 1-dimensional chain. We give the following result: 

Theorem 2 The 2-local Hamiltonian on a 1-dimensional chain of 13-state particles is QMA- 
complete. This in turn implies that the 8-local Hamiltonian on a 1-dimensional chain is QMA- 
complete. 

Certainly an intriguing question left open here is whether the locality can be reduced and, in 
particular, whether the 2-local Hamiltonian on a 1-dimensional chain is QMA-complete. 

Similar results to these have been independently obtained by Aharonov, Gottesman and Kempe 

[?]■ 



2 



2 Preliminaries 



2.1 The Model of Adiabatic Computation 

The Adiabatic Theorem, stated here, is the foundation of adiabatic computation. 

Theorem 3 The Adiabatic Theorem (as adapted from [?] and quoted from [?] ) Let 

Hi n it and Hu na i be two Hamiltonians acting on a quantum system and consider the time- dependent 
Hamiltonian H(s) = (1 — s)Hi n i t + sHf ina i. Assume that for all s, H(s) has a unique ground state. 
Then for any fixed 5 > 0, if 

rp ^ q ( \\H final ~ Hj n jt\\ 1+S \ , . 

- U \eSmin sem {A^(H(s))})> W 

the final state of an adiabatic evolution according to H for time T (with an appropriate setting of 
global phase) is e-close in l2-norm to the ground state of Hfi na i. The matrix norm is the spectral 
norm \\H\\ = m&x w \\Hw\\/\\w\\. 

6 will be a contstant and the constant in the 0, will go to infinity as 5 goes to 0. We use the model 
of adiabatic computation as described in [?]: 

Definition 1 A k-local adiabatic computation AC(n,d, Hi n n, Hfi na i,e) is specified by two k-local 
Hamiltonians, Hi n a and Hfi na i acting on n d-dimensional particles, such that both Hamiltonians 
have unique ground states. The ground state of Hi n n is a tensor product state. The output is a 
state that is e-close in l2-norm to the ground state of Hji na [. Let T be the smallest time such that 
the final state of an adiabatic evolution according to H(s) = (1 — s)Hi n a + sHfi na i for time T is 
e-close in l s -norm to the ground state of Hfi na i. The running time of the adiabatic algorithm is 
defined to be T ■ max s ||iT(s)||. 

Thus, the running time of the computation will be upper bounded by the norm of the Hamiltonians 
and the lower bound given in Inequality [TJ 

It is possible in some cases to restrict our attention to a subspace S. Suppose that H(s) leaves S 
invariant. H$(s) is the restriction of H(s) to S. If we start the adiabatic evolution at a state inside 
S then an adiabatic evolution according to H$(s) is identical to an adiabatic evolution according 
to H(s). Therefore, we can use AHs(s) in Theorem [3] instead of AH(s). 

2.2 The Class QMA and the Local Hamiltonian Problem 

The class QMA is defined in terms of promise problems defined by a pair (L yes ,L no ) of disjoint 
sets of strings corresponding to YES and NO instances of the problem. The input is guaranteed 
to be in L yes U L no and solving the problem entails determining whether a given input string x is 
in L no or L yes . B is defined to be the Hilbert space of a single quibit. 

Definition 2 (QMA) Fix e = e(\x\) such that e = 2~ n (\ x \\ Then, a promise problem L G QMA 
if there exists a quantum polynomial time verifier V and a polynomial p such that 

V* G L yes , 3|0 G B®rtW\Pr(V(\x), |0) = 1) > 1 - e 

V* G L no ,V|0 G B^\Pr(V(\x), |0) = 1) < e 
where Pr(V(\x), = 1) denotes the probability that V outputs 1 given \x) and 
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The first known QMA-complete problem is the local Hamiltonian problem and is a natural 
analog of SAT . 

Definition 3 We say that an operator H : B® n — > B® n on n qubits is k-local if H can be expressed 
as the sum of terms, where each term is a Hermitian operator acting on at most k bits. 

Definition 4 (The (promise) problem k-local Hamiltonian ) A k-local Hamiltonian on n 
qubits H = Y<j =1 Hj with r = poly(n). Each Hj has a bounded operator norm \\Hj\\ < poly(n) and 
its entries are specified by poly{n) bits. In addition, we are given two constants a and b with a < b. 
In YES instances, the smallest eigenvalue of H is at most a. In NO instances, it is larger than 
b. 

2.3 Previous Results on Spectral Gaps 

We make use of some previous results on spectral gaps which we state here in a form that is 
particularly suited for our purposes. Define P r to be the r x r matrix of the following form: 

/ \ -| ••• 0\ 

-| i -i o ••• 

0-1 1 -3 '•• 

V o ••• " o -I \ I 

Lemma 4 Consider a Hamiltonian H = Hi + H2 acting on a subspace S. Consider a basis for 
S, 71, ... ,7,.. Suppose that when H\ is restricted to S and expressed in this basis, the matrix is 
P r . Suppose also that when H2 is restricted to S and expressed in this basis, the matrix is diagonal 
with non-negative integer entries, at least one of which is nonzero. Then the lowest eigenvalue of 
H when restricted to S is f2(l/r 4 ). 

For the proof of this lemma, we make use of the following lemma proven by Kitaev (Lemma 
14.4 in [?]). 

Lemma 5 Let H\ and H2 be two Hamiltonians with ground energies a± and 02, respectively. Sup- 
pose that for both Hamiltonians the difference between the energy of the (possibly degenerate) ground 
space and the next highest eigenvalue is larger than A, and that the angle between the two ground 
spaces is 9. Then the ground energy of H\ + H2 is at least a\ + 02 + 2Asin 2 (#/2). 

Proof of Lemma [4j Suppose that the diagonal entries are all non-zero. Since the entries are 
integral, the lowest eigenvalue of H2 is at least 1. Since H\ and H2 are both positive semi-definite, 
the ground energy of H\ + H2 is at least 1. Now suppose that H2 has at least one non-zero entry 
It's spectral gap is at least 1 since it has at least one non-zero entry on its diagonal. Using standard 
techniques, one can show that the spectral gap of H\ is at least f2(l/r 2 ). The ground energies of 
H\ and H2 are both 0. Furthermore since the ground state of H\ is a uniform superposition of all 
the basis vectors and H2 has at least one non-zero entry, the cosine of the angle between the two 
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ground spaces is at most 1 — 1/r. Invoking Lemma [5l we have that the ground energy of Hi + H2 
is at least J7(l/r 3 ). ■ 

The next lemma is instrumental in providing a spectral gap for the adiabatic computation: 

Lemma 6 Let Hf ina i be P r for some r. Let Hi n a be an r x r matrix of the following form: 

/ ••• \ 
1 ••• 

V° ••• !/ 

Define H(s) = (1 — s)Hi n n + sHfi na [. Then the spectral gap of H(s) is S7(r~ 2 ) for all s £ [0, 1]. 
This lemma is proven in [?]. 

3 A One-dimensional Universal Quantum Adiabatic Computer 

The proof of the following theorem shows that we can simulate a quantum circuit that uses L gates 
on n qubits with an adiabatic computation on a chain of (L + 2) ten-state particles on a line. 

Theorem 7 Given a quantum circuit on n qubits with L two-qubit gates implementing a unitary 
U, and e > 0, there is a 2-local adiabatic computation between nearest-neighbor particles on a line 
AC(L + 2,W,Hi n it,Hfi na i,e) whose running time is polynomial in L and 1/e and whose output is e 
close to U\0 n ). Moreover, Hi n n and LLjinai can be computed by a polynomial time Turing machine. 

We first show a consturction that uses particles with thirteen states. We will add one extra 
state at the very end and then reduce the number of states by four by identifying pairs of states. 
The final construction will use ten states. Our proof makes use of the observation in Section 12.11 
that one need only prove a lower bound on the spectral gap of H(s) when restricted to a subspace 
as long as the adiabatic computation starts in that same subspace. 

We follow the convention in [?] which assumes without loss of generality that the quantum 
circuit to be simulated has a particular layout of the gates. It consists of a sequence of R rounds. 
Each round is composed of n nearest neighbor gates. The first gate in each round is a one-qubit 
gate applied to the first qubit. For i = 2, . . . , n, the i th gate is a two-qubit gate applied to qubits 
i — 1 and i. It will be convenient to assume that the first round of gates consists entirely of identity 
gates. This will just serve to increase the number of rounds by one. The value of R is adjusted 
accordingly. The construction in [?] also has a sequence of n indentity gates at the end of each 
round but we omit that here. The total number of gates then is L = nR. They observe that any 
circuit can be transformed to fit this format by adding additional identity and swap gates. 

We now define a Hamiltonian H prop which will enforce the propogation of the quantum com- 
putation. In this section Hn na i will simply be H prop . The ground state of the Hu na i will be a 
uniform superposition of some number of different configurations which we call templates. Each 
template represents a subspace of 2 n possible states and the particular state within a given template 
will encode the state of the computation at some particular point in time. We will think of these 
templates as a sequence and that the state of the system changes through time from one template 
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to the next. In fact, the ground state of Hjmai will be a superposition of these different snapshots 
of the computation. We start by giving an overview of H prop . In the course of this description, we 
will describe the states for each particle and their significance. 

We have a chain of L + 2 particles, where L is the number of gates in the quantum circuit. The 
particles will be labelled through L + 1. We will think of the particles as organized into L/n 
contiguous blocks of n particles with a single particle on either end (particles and L + 1). The 
first block of particles are 1 through n, the second block is n + 1 through 2n, etc. At any point in 
time, there are n particles whose state represent the state of the computation. We will call these 
the computation particles. These are either all contiguous or possibly separated by a single control 
particle somewhere in the middle. In general, H prop will enforce that the computation quibits shift 
from the left to the right. When the computation particles are all aligned within the boundaries 
of a single block (i.e. are located in positions in + 1, ...,(« + l)n for some i), then the gates are 
applied for the i th round. 

There are 13 possible states for the qubits. 

• F: F stands for finished. Particles in this state are to the left of the computation particles 
and will no longer change state. 

• N: N stands for new. Particles in this state are to the right of the computation particles and 
have not yet been reached. 

• L: L stands for left. A particle in this state is used for control and indicates that the particle 
in the control state is propogating to the left. 

• Tr, Tl : Tl and Tr stand for turn left and right, respectively. A particle in this state is used 
for control and indicates that the control particle will start moving in the opposite direction. 
Tl indicates that it will start moving to the left and Tr indicates that it will start moving to 
the right. 

• Qo 5 Qi>Bo,Bi : These states represent qubits. Particles in these states are computation 
qubits. The subscript represents the state of a qubit in the computation. The B states are 
to the left of the particle in the control state and the Q states are to the right of the particle 
in the control state. 

• Ro,Ri : R stands for right. A particle in this state is a computation particle as well as a 
control particle. The subscript represents the state of a qubit in the computation. A particle 
in this state also indicates that the control particle is moving to the right. 

• Go,Gi : G stands for gate. A particle in this state is a computation particle as well as a 
control particle. The subscript represents the state of a qubit in the computation. A particle 
in this state also indicates that the control particle is moving to the right and a gate is being 
applied. 

We will explain the basic construction of H prop by showing an iteration in which the control 
particle propogates to the right and then to the left. At the end of the iteration, the computation 
particles will have shifted one position to the right. At each step, we indicate the component of 
H prop that enforces the propogation. If A, B, X, Y are possible particle states, we will define 

H'ab-xy = \AB){AB\ hl+1 + \XY){XY\ hl+1 - \AB)(XY\ iji+1 - \XY)(AB\ hl+1 . 
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The superscript i indicates that the Hamiltonian is applied to particles i and i + This particular 
Hamiltonian enforces the condition that if a state is an eigenstate with eigenvalue 0, then for each 
configuration with AB in positions i and i + 1 there must be a state of equivalent amplitude with 
AB replaced by XY. 

It will be useful in our exposition for each H l AB ^ XY to define two companion tranformations. 
These correspond to transforming any state with AB in locations i and i + 1 by changing the AB 
to XY. This is represented by J 7 , where the T stands for propogation in the forward direction. 
Similarly B will correspond to replacing XY by AB which represents propogation in the backward 
direction. 



AB^XY 



\XY){AB\ iti+1 , B 



AB*—XY 



\AB)(XY\i 



i+1- 



These will not be part of the final Hamiltonian but will be helpful for describing certain states. 

For the time being, we will not use that gate states (Go and Gi). This means that the state 
of the represented quibits will not change and there is no change of state between {Qq,Rq} and 
{Qi,Ri}. For ease of notation, we will omit the subscripts. For example, Hbl^lq actually 
represents Hb l^lq + Hb\L^lQ\- Similarly the expression Hrq^br would be the sum of four 
different terms for the possible subscripts on the pair of states. We will use the term template to 
be a representation of a set of states that all have the same letter value with the subscripts on the 
computation bits unspecified. 

We start with the particles in the following template: 





F 




F 


Tr 


Q 


Q 


Q 


Q 


Q 


N 


N 




N 




where n = 5. Let particle i b 
template: 


e the location of the Tr state. Hj, 


rQ<-»- 


FR will result in 




F 




F 


F 


R 


Q 


Q 


Q 


Q 


N 


N 




N 




For j = i + 1 through i + n — 


1 we have the Hamlitonian H 3 R q^ 


BR which results in 




F 




F 


F 


B 


B 


B 


B 


R 


N 


N 




N 




Then H 1 r + £^ BTl , results 


in 




























F 




F 


F 


B 


B 


B 


B 


B 


T L 


N 




N 





The presence of the particle in the Tl (turn left) state indicates that propogation will switch from 
moving right to moving left. Ht^^ ln , results in 





F 




F 


F 


B 


B 


B 


B 


B 


L 


N 




N 




The presence of the particle in state L indicates propo§ 
down to i + 1 until 


Ration to the left via H 3 bl ^ L q for j 




F 




F 


F 


L 


Q 


Q 


Q 


Q 


Q 


N 




N 





n 



Finally, H 1 fl<->ft r brings us back to the original template with the computation particles shifted 
one position to the right: 



F 




F 


F 


T 


Q 


Q 


Q 


Q 


Q 


N 




N 
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There are a total of 2n + 3 configurations in one iteration that shifts the states one position to the 
right. 

The application of the gates is triggered when the computation qubits align within the bound- 
aries of a set because we have a slightly different Hamiltonian at these locations. At the locations 
i such that i is a multiple of n, we have H Fr q^ fg instead of H Fr q^ fr . In addition, we apply 
the qubit gate to Qo/i- If the Q is in position i, then this is the i th gate in the computation and 
we denote the 2x2 unitary matrix associated with this 1-bit gate as C/j. We will express this 
Hamlitonian matrix over the subspace spanned by {TrQo,TrQi, FGq, FG\] 



H. 



T R Q^FG 



I -Ui 
-Uj I 



FtrQ^fg 



Ui 




B T B Q*-FG 





uj 








Similarly for the 2-qubit gates, we apply the i th gate to GQ in locations i — 1 and i and the 
result is a state of the form BG. We express this Hamiltonian as a matrix over the basis 

{GqQq, GqQi,GiQo, GiQi, BqGq, BqG\, B\Gq, B\Gi}, 



where the i th gate is expressed by the 4x4 unitary martix Uf 



H, 



GQ^BG 



I 

-u 



-Ui 

I 



<ri—l 

^GQ^BG 



Ui 




B. 



:-l 

GQ^BG 



U 








Finally we need the term that turns the propogation of the control state when G reaches the 
rightmost end of the computation bits. For i mod n = 0, we have Hqn^bt l instead of Hrn^bt l - 
This just preserves the value of the represented qubit as before. To summarize, 



H, 



prop 



prop' 



For i mod n/0 



H 



prop 



l RN^BT L 



+ H bl ^ L q + Hq Q ^ bg + H fl ^ FTr 



T L N^LN 



For i mod n = 0, we replace H 1 TrQ ^ fr with W TrQ ^ fg and H rn ^ BTl with H 1 gn ^ BTl , to get 



H 



prop 



= H. 



T R Q~FG 



RQ^BR 



GN^BT L 



T L N^LN 



+ Hbl^lq + H G q^ bg + H fl ^ FTr 



Finally, if i = L, we remove H gn ^ BTl to end the sequence. We can define T by using the corre- 
sponding T term in place each H prop term. Similarly we can define B by using the corresponding 
B term in place each Hprop term. 

Now we would like to define the target ground state for Hfi na i. In order to define the target 
ground state, we define a series of templates. These templates are linearly ordered. The label t for 
each template ranges from through T = (2n + 3)(L — n) + n + 1. The template is represented 
by a string of length L + 2 indicating the state for each particle, omitting the subscript for any 
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Condition 


Template 


Forward 


Backward 








(i+i) 


••(m+i) 




(n+i+l)- 
























if (not 


boundary) 






j = 


F • 


• • F Tji 


Q • 


• Q 




AW-- 


■N 


TrQ 


-> FR 


FL <- 


-FT R 


















if (boundary) 






















TrQ 


-» FG 






i mod j = 1 


F • 


■■F F 


RQ--- 


QQ-- 


Q 


NN- ■ 


■N 


RQ 


-> BR 




<- Fi? 


i mod n = 0, j = 1 


F • 


■■F F 


GQ--- 


QQ-- 


Q 


N N- 


■N 


GQ 


-> 




«- FG 


i mod n^O, Kj<ra 


F • 


■■F F 


B ■■■BR Q-- 


_Q 


N N- 


■N 


RQ 


-» 5F 


RQ- 


- Si? 








j 


m-j 


-i 














j mod n=0, Kj<n 


F- 


■■F F 


B ■ ■ ■ BG Q-- 


_Q 


N N- 


■N 


GQ 


-» 


GQ 


- 








3 


m—'j 


-i 














i mod n^O, j = n 


F- 


■■F F 


B- 


■■BR 




N N- 


■N 


RN 


-> ST 


RQ- 


- BR 


i mod n = 0, j = n 


F- 


■■F F 


B- 


■BG 




N N- 


■N 


GN 


-> 


GQ- 


- 






















if (not 


boundary) 


j = n + 1 


F- 


■■F F 


B 


■■■B 




T L N- 


■N 


T L N 


LiV 


RN < 


-BT L 






















i/ (boundary) 






















GN * 


-BT L 


j = n + 2 


F- 


■■F F 


B 


■ B 




LN ■ ■ 


■N 


BL 




T L N 


«- LiV 


n + 2<j<2n + 2 


F- 


■■F F 


B---BL Q-- 


■_Q 






BL 




BL 


-LQ 








m-j 


j- 


i 


L-i 


— m 










j = 2n + 2 


F- 


■■F F 




■Q 




QN ■■ 


■N 


FL - 


^FT R 


LQ- 


- BL 








m—l 

















Figure 1: Templates 



computation bits. In this section, the templates will always have n computation bits. Thus, 
the remaining degrees of freedom represented by the omitted subscripts results in a subspace of 
dimension 2 n . We say that a term Hab^xy m H prop applies in the forward direction to a template 
if AB appears in string corresponding to the template. This means that J-ab— >xy has a non-zero 
result when applied to a state in the corresponding subspace for that template. The result will 
be a state in the template with AB replaced by XY. Similarly Hab^xy applies in the backward 
direction if XY appears in the template. Bab^xy applied to a state in that template results 
in a state in the template with XY replaced by AB. Figure [3] shows a chart of the templates. 
The string representing the templates will depend on the value of t as indicated by the column 
labelled Conditions. We also include a column that indicates for each template the set of terms in 
Hp r0 p that apply in the forward direction and a column for the terms that apply in the backwards 
direction. The condition Boundary holds for a given term if the pair straddles a boundary between 
blocks of n particles. That is, if it applies to particles p and p + 1 and p mod n = 0. For each t, let 
i = [t/(2n + 3)J and j = t mod (2n + 3). 

Lemma 8 For each t G {0, . . . ,T — 1}, exactly one term in H prop applies to the t th template. 
Furthermore, when Tp TO p is applied to a state in that template, it results in a state in the next 
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template in the sequence. For each t S {1,...,T} ; exactly one term in H prop applies to the t th 
template in the backward direction. Furthermore, when B prop is applied to a state in that template, 
it results in a state in the previous template in the sequence. B is zero on a state in template and 
T is on a state in template T. 

Proof: The proof consists of verifying that for each template in Figure El exactly one term applies 
in the forward direction (and in only one place) and exactly one term applies in the backward 
direction (again in only one place). When the term is applied in the forward direction, the result 
is the next template in the sequence and when the term is applied in the backward direction, the 
result is the previous template in the sequence. Template is T R Q n N L - n+1 and there are no 
terms in H prop that apply in the backward direction. The T th template is F L ~ n+1 B n ~ 1 GN . Since 
^gn*-*bt * s removed from H prop , there is no term in H prop that applies in the forward direction. 
■ 

Each template with n computation particles corresponds to a particular point in the compu- 
tation to be simulated. However, the correspondence is not one to one. A particular point in the 
computation can be represented by many templates. The t th template corresponds to the point in 
the computation after g(t) gates have been performed, where g(t) is the number of templates in the 
sequence from through t containing a particle in a G state. As each new template with a particle 
in a G state is reached, one more gate is performed. If i = \t/(2n + 3) J and j = t mod (2n + 3), 
then 

i + min{n, j} if i mod n = 0, 
n\i/n\ if i mod n ^ 



g(t) 



Define \4>t) to be the state in the t th template such that the superposition of states defined by 
the 0/1 subscripts of the computation bits corresponds to the superposition of states in the circuit 
after g(t) gates have been applied, assuming that the input to the quantum computation is |0). 
Thus, |0o> = T(Q ) n N L - n+1 . 

For each t € {0, . . . , T — 1}, T prop \4>t] = \<f>t+l) and ^(fir) = 0- Also, for each t £ {1, . . . , T}, 
B\tpt) = l^t-i) an d B\(f>o) = 0. Now define S to be the subspace spanned by the states \4>t) for all 
t E {0, . . . , T}. These states form an orthonormal basis of S. S is closed under Hprop. The matrix 
representation of H prop in the \4>) basis is exactly Pt+i- Recall the definition of P r from Section 
12.31 The unique ground state of Hf ina [ ( = H prop ) is 

1 T 

W final) = -==$^|0t). 
V 1 + 1 i=0 

Now we need to define Hi n n\ 

Hinit = {I — \Tr)(Tr\)q. 

Note however, that while Hi n n in its current form can distinguish between |0q) and the other |0)'s, 
it does not ensure that the initial configuration starts out in S. To address this problem, we add 
an additional state S (for Start) and some extra constraints. The initial configuration will be 
TrS 11 N L ~ n+1 . We will add a constraint to Hi n u that enforces the condition that if particle is in 
state Tr, then particle 1 is in state S. This is achieved by adding in \TrX){TrX\o, where we sum 
over all states X such that 1/5. Similarly for i = 1 through n — 1, we add the constraint that 
if particle i is in state S, then i + 1 is also in state S. We also add the constraint that if particle 
n is in state S, then particle n + 1 is in state N. Finally for all i = n + 1 through L, we add the 
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constraint that if particle i is in state N, then particle i + 1 is in state N. All of these constraints 
are added into i?i n tf- 

We need to alter H prop slightly in order to work with the S states. In H prop , we replace Hj, r q<^ fg 
with Hj, r s^,fq - (Recall that the input is all O's). Then for i = 1 through n— 1, we replace H G q^ bg 
with Hq S ^ BqG . Templates 1 through n are now of the form F{B y G S n - j - l N L - n+1 . Lemmal 
still holds. Furthermore, ensures that there is only one eigenstate with eigenvalue and it is 
our desired initial configuration. 

Lemma 9 The matrix representation of Hinit restricted to S and expressed in the basis of \4>t)'s 
is a diagonal matrix. The first diagonal entry is and the others are 1. 

Proof: All templates start with F except the first one which starts with T. All the templates satisfy 
the other conditions in Hi n it- ■ 

To establish the spectral gap of Hs (s), we simply invoke Lemma[6]which says that the spectral 
gap of the restriction of H(s) to Sq satisfies A(i?s(s)) = f2((Lra)~ 2 ). 

Using the Adiabatic Theorem, we get that the running time of the algorithm will be 0(\ \H(s)\ \e~ s (Ln] 
Note that ||iJ(s)|| is O(L) which gives an overall running time of 0(e _l5 L(Ln) 4+2 ' 5 ). However, this 
only produces a final state that is e close to \(f> final)- This is a superposition of the |^>t)'s and 
only \(j>r) encodes the desired output. This can be corrected by adding L/e identity gates to the 
computation so that only a fraction of e of the j)'s will encode partial points in the computation. 
This makes the final running time 0(e~( 5+3<5 )-L(Ln) 4+25 ). 

In order to reduce the number of states from 14 to 10, we observe that we can identify pairs 
(N,F), (T R ,T L ), (B ,Qo) and and Lemma [8] still holds. 



4+2<5> 



4 8-local Hamiltonian on a Chain is QMA-complete 

We now turn to the problem of the promise local Hamiltonian. We will continue to work with 
13-state particles as in the previous section. (We will not use the 5 state added at the end of the 
previous section). The ground state will now encode the computation performed by a quantum 
verifier V which works on input \x£) for input x and witness £. The total length of the input is n 
qubits. We will construct a Hamiltonian such that if there exists a witness £ such that V accepts 
with probability at least 1 — e on input |0£), then the lowest eigenvalue of H will be less than e. 
However, if for every £, V accepts with probability at most e, then the lowest eigenvalue will be at 
least 1 over a polynomial in n and L. 

For this problem, we need to show that there is a large spectral gap over the entire space of the 
particles, not just when restricted to a particular subspace. The Hamiltonian H will consist of five 
components: 

H — Hprop "I" Hinit H ou i -\- H va nd -\- H\ e g a \. 

Hprop w iH the same as in the previous section (with the stipulation that we do not use the changes 
made at the end of the previous section to incorporate the additional state S). 

We will use the term template here to refer to any L + 2-character string over the alphabet 

{F, N, B, Q, L, R, G, T R , T L }. 
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Figure 2: Graph indicating set of valid templates. 

As before, each template will designate a subspace according to the unspecified 0/1 subscripts for 
the computation states. This subspace will have dimension 2 m if the template has m computation 
bits. 

We will use various hamiltonians to enforce that only certain templates will be valid or legal. 
We start with a set of consraints enforced in H va ud- We will enforce through H va nd that a template 
must have a form which is a string of length L + 2 and is specified by one of the three regular 
expressions: 

F + B*(R + G + L)Q*N + , F*T L Q + N + , F + B + T R N. 

Another way to denote these constraints is that a string corresponding to a valid template must be 
a path in the graph in Figure d] from the Start node to the End node consisting of L + 2 internal 
nodes. R, L, G are grouped together for clarity. A path through this node can use either R, L or 
G. 

These constaints are enforced by having a set of allowed pairs where an edge in the graph 
corresponds to an allowed pair. We will need some additional constraints to enforce that the first 
character must be F or T R and that the last characers must be N. For ease of notation, we will 
omit the subscripts. Therefore, the pair RQ represents four possible pairs for the four possible 
combinations of subscripts. When summing over a set, summing over all possible combinations of 
subscripts will be assumed. 

S = {FF, FB, BB, FR, FL, FG, BR, BL, BG, RQ, GQ, LQ, 
QQ, QN, NN, RN, GN, LN, FT R , T R Q, BT L ,T L N} 
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Another way to express the set S is that it is the set below, where the characer X can be either R, 
L, or G: 

S = {FF,FB,BB,FX,BX,XQ,QQ,QN,NN,XN,FT R ,T R Q,BT L ,T L N} 

There is a one-to-one correspondence between the above pairs and the edges in Figure [H 
For any i from to T — 1, we have 

H Llid = I - \ a P)( a P\i,i+i- 

a/3eS 

We then sum these together and add some additional constraints at the beginning and end of the 
chain. 

L 

H valid = (I-\F)(F\ - \T R ){T R \) + (I-\N)(N\) L+1 + J2Hi egal . 

i=0 

Note that H va ud is a diagonal matrix with non-negative integers along the diagonal. 

Lemma 10 Consider a state \<j>) contained in a subspace specified by a template. H va nd\4>) = X\4>) 
for some non-negative integer A. If A = 0, then the template must have a form specified by one of 
the following three regular expressions: 

F + B*(R + G + L)Q*N + , F*T L Q + N + , F + B + T R N. 

Proof: If H va iid\4>) = 0, then every consecutive pair of characters in the the template containing 
\<ft) must be an allowable pair. Furthermore, the first character must be F or T R and the last 
character must be N . The graph in Figure H] has every allowable pair labelled as a directed edge. 
Furthermore, the graph enforces that any path from Start to End must begin with an F or a T R and 
must end with an N. Therefore, a valid template must correspond to a path that starts at the Start 
node, ends at the End node and traverses L + 2 intermediate nodes. The length of the template is 
enforced by the physical length of the chain of particles. The set of all such paths correspond to 
the three regular expressions above with the constraint that the length must be L + 2. ■ 

We will label each valid template with a pair (m, t). m will designate the number of computation 
bits in a template. We describe here how to determine t for a particular m. We describe how to 
determine two inetegers i and j and then let t = i(2m + 3) +j. The labelling will have four distinct 
cases, depending on the form of the template: 



F k Ti y Q m N L ~ m ~ k+1 


i <— k 


3^0 


F k+1 B L (G + R)Q m - l ~ 1 N L - m - k+l 


i <— k 




pk+l j^mrp^jy L—m— k+1 


i <— k 


j <— m + 1 


pk+1 j^m—l JjQl jyL—m—k 


i <— k 


j <- m + 2 + / 



Conversely, Figure H] shows for a given pair (m, t) the form of the templates corresponding to 
that pair, where i = |_t/(2m + 3)J and j = t mod (2m + 3). Note that for each pair, there is exactly 
one template unless 1 < j < m in which case there are two templates, depending on whether the 
particle in the control state is in a G or an R state. Define T m = (2m + 3)(L — m) + m + 1. Also 
note that for a given m, and any t E {0, . . . , T m }, (m, t) corresponds to a valid template, but having 
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Condition 


Template 


Forward 


Backward 






0---i 


(i+i) 


■■(m+i) 


(m+i+l)- 


(L+l) 




















if (not 


boundary) 






3 = 


F • 


■■FT R 


Q ■ 


■ Q 


NN-- 


■N 


TrQ 


-> FR 




-FT R 
















if (boundary) 




















TrQ 


-» FG 
























i/ (not 


boundary) 


3 = 1 


F ■ 


■■F F 


RQ--- 


QQ---Q 


NN ■ ■ 


■N 


RQ 


-» BR 


TrQ 


*- FR 




















if (boundary) 


3 = 1 


F ■ 


■■F F 


GQ--- 


QQ---Q 


NN ■ ■ 


■N 


RQ 


-» BG 


TrQ 


^FG 


1 < j < n 


F ■ 


■■F F 


B ■ ■ ■ B R Q---Q 


N N ■ 


■N 


RQ 


-»■ BR 


RQ 


- BR 








j 


m—j—l 














1 < j < n 


F ■ 


■■F F 


B ■ ■ ■ BG Q---Q 


N N ■ 


■N 


GQ 


-> BG 


GQ 


- BG 








j 


m-j -I 




























if (not 


boundary) 






j = n 


F ■ 


■■F F 


B- 


■BR 


N N • 


■N 


RN 


-» BT 


RQ 


- BR 
















if (boundary) 






j = n 


F ■ 


■■F F 


B- 


■BG 


N N ■ 


■N 


GN 


-> BT 


GQ 


- BG 




















if (not 


boundary) 


j = n + 1 


F ■ 


■■F F 


B 


■■■B 


T L N-- 


■N 


T L N 


LN 


RN i 


-BT L 




















if (boundary) 




















GN * 


-BT L 


j = n + 2 


F ■ 


■■F F 


B 


■ B 


LN •■ 


■N 


BL 


-> LQ 


T L N 


«- LN 


n + 2<j<2n + 2 


F ■ 


■■F F 


B ■ ■ ■ B L Q -Q 




__/V 


BL 


-> 


BL 


-LQ 








m-j 


3-1 


L-i 


— m 










j = 2n + 2 


F ■ 


■■F F 




■Q 


QN- 


■N 


FL - 




LQ 


- BL 








m—l 















Figure 3: Generalized Templates 



t > T m implies that the last particle will not be in state N, which makes the template invalid. If 
a template is labelled (m, 0), then we will call it an initial template. If it is labelled (m,T m ), then 
we will call it a final template. 

Figure S] also shows for each template which terms in H prop apply in the forward direction 
and which apply in the backward direction. The condition (boundary) indicates whether a pair 
straddles two blocks of particles. That is, if a pair of particles are located in positions i and i+1 
and i mod n = 0. 

Lemma 11 Consider a valid template T labelled with (m,t). At most one term in H prop applies to 
T in the forward direction. Furthermore, when T is applied to a state in that template, the result 
is or a state in an (m,t + \)-template. At most one term in H prop applies to T in the backward 
direction. Furthermore, when B is applied to a state in that template, the result is or a state in 
an (m,t — 1) -template. 
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Proof: The proof consists of verifying that for each template in Figure HI at most one term applies 
in the forward direction (and in only one place) and exactly one term applies in the backward 
direction (again in only one place). When the term is applied in the forward direction, the result 
is the next template in the sequence and when the term is applied in the backward direction, the 
result is the previous template in the sequence. ■ 

Lemma 12 Consider a state \<f>) contained in the subspace corresponding to a template T. IfTis 
an initial template, then B\<f>) = 0. If T is a final template, then J-\4>) = 0. 

Proof: There is no term in H prop that applies to an inital template in the backward direction. 
Furthermore, since Hq N ^ bt is removed from Hp rop , there is no term in H prop that applies to a 
final template in the forward direction. ■ 

We define another Hamiltonion H[ ega i which penalizes any template for which there is no term 
in H prop that applies in the backward direction or for which there is no term that applies in the 
forward direction (unless it is an inital or final template, respectively). Using the table in Figure [4j 
we want to forbid pairs RN and FR from crossing a boundary. We will also forbid pairs GQ and 
BG from crossing a boundary. In addition, we want to forbid pairs GN and FG unless they cross 
a boundary. 

For any i such that i mod n = 0, we have 

Hi ega i = \RN)(RN\i + \FR)(FR\i + \GQ)(GQ\ t + \BG)(BG\ t . 

For any i such that i mod n ^ 0, we have 

H} egal = \GN){GN\ l + \FG){FG\ t . 

As usual, summing over all the subscripts of the computation states is assumed. Finally, we sum 
these together 

L 

Hlegal = H\ egal . 

We say that a template is legal if any state |0) in that template has Hi ega i\(p) = 0. Otherwise, 
Hi ega i\(j)) > 1 and we say the template is illegal. 

Lemma 13 Consider a template T that is both valid and legal. If T is not a final template, then 
there is a term in H prop that applies to T in the forward direction. Similarly, if T is not an initial 
template, there is a term in H prop that applies to T in the backward direction. 

Proof: The proof consists of the obervation that any template in Figure 2] for which there is no 
term in H prop that applies in the forward direction is made illegal by Hi ega i. Similarly, any template 
for which there is no term in H prop that applies in the backward direction is made illegal by Hi ega i ■ 
■ 

We can think of all the valid templates as nodes in a directed graph. There is a directed edge 
from templates T to T if applying T to some state in T results in a state in T . (By definition 
then, applying B to a state in T results in a state in T). We will call this graph the template graph 
and will refer to nodes and templates interchangeably. Lemmas [TT1 [T2l and [13] imply that this graph 
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consists of a set of disjoint chains. All the nodes in a chain correspond to templates with the same 
number of computation particles. Furthermore, the starting node in a maximal chain is either an 
initial template or an illegal one. The last node in a maximal chain is either a final template or an 
illegal one. 

Lemma 14 There is exactly one chain in the template graph that contains no illegal nodes. Fur- 
thermore, templates in this chain have n computation particles. 

Proof: Consider a maximal chain with no illegal nodes. This chain must begin with an initial 
template and end with a final one. The initial template has the form TrQ" 1 N L ~ m+1 for some m. 
By the forward rule T R Q -> FG, the next node in the chain is FGQ m ' 1 N L ' rn+1 . If m < n, then 
m — 1 applications for the forward rule GQ — > BG will result in FB m ~ 1 GN L ~ m ~ 1 . Since m < n, 
this will result in a pair GN in locations m and m + 1 which does not straddle a block boundary. 
This is made illegal in Hi ega i. If m > n, then n — 1 applications of the forward rule GQ — > BG to 
FGQ m ~ 1 iV L - m+1 will result in FB^GQ™' 11 ^-™* 1 . This will result in a pair GQ in locations 
n and n + 1 which is also disallowed in Hi ega i. 

Finally, if m = n, there is exactly one chain from an initial template to a final one (because 
there is exactly one initial template). The application of the forward rules do not result in any 
illegal tempaltes. ■ 

It will be convenient at this point to define the remaining two terms in H. The input to the 
quantum verifier will be n qubits. n\ qubits will be ancillary qubits that are initialized to and n 2 
qubits will be the witness £. n\ + n 2 = n. We force x = with the following Hamiltonian: 

m 

Hinput = X! \Qi)(Qi\i + + \Gl)(Gi\i + \Bi)(Bi\i. 

i=i 

We will assume that the output will be present in the rightmost qubit in of the computation. 
Therefore, we will have H out defined to be \Go)(Gq\l- Observe that H input and H out are both 
closed over the subspace defined by each template. 

A maximal chain in the template graph defines a subspace which is just the subspace spanned 
by all the subspaces defined by the templates along the chain. H prop is closed over the subspace 
defined by any maximal chain in the template graph. All the other terms in H are closed over the 
subspace defined by each template. 

Define Si ega i to be the subspace defined by the unique chain containing only legal nodes. Let 
^legal De ^ ne orthogonal subspace to Si ega i. Since H is closed under Si ega i, it is also closed under 
^legal ana - an y eigenvector of H must be contained in Si ega i or S^ gal . 

Lemma 15 Any eigenvector of H in S^ gal will have an eigenvalue of at least Q(1/L A ). 

Proof: Define S va ud to be the subspace spanned by all the valid templates. Since H is closed over 
Svalid, h is also closed under the orthogonal space S^ aUd . Any state eigenstate of B vo xid i n S^ alid 
will have an eigenvalue of at least 1. Since the remaining terms in H are positive semi-defininte, 
any eigenstate of H in S^nd will also have an eigenvalue of at least 1. 

Now we can carve up H va i id into the subspaces defined by the maximal chains and H is closed 
on each such subspace. We focus on an arbitrary such maximal chain (except the one containing 
only legal nodes) and the subspace it defines. The chain goes from an (m,ti) template to an 
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(m, ^-template for some m, t± and t-i- To specify a state in the first template, we specify an 
m-bit x string which will determine the subscripts of particles in the computation states. We call 
this state \(f>x,m,ti) ■ The set of these states for all x forms a basis of the first template. Since T 
is unitary, when we apply T to all the |0}'s, we get a basis of the next template in the chain. 
Applying T i — 1 times gives a basis of the i th template in the chain. We will focus on a sequence 
), where .P|<Az,m,ti) = \<j>x,m,t±+i) ' Tne subspace spanned by these states is 
closed under H prop + Hiegai- Furthermore H prop when restricted to this subspace and expressed 
in the basis of 0's is P r , for r = t% — t\ + 1. Hi ega i when expressed in this basis is diagonal with 
non-negative integer entries. Furthermore, we know that there is at least one positive entry on the 
diagonal because the chain has at least one illegal node. By Lemma [U we know that the lowest 
eigenvalue of any eigenstate in the subspace defined by the chain must have an eigenvalue that is at 
least r2(l/(T n ) 3 ). Since the remaining terms in H are all positive semi-defininte, the lower bound 
holds for H as well. ■ 

We will prove the following theorem: 

Theorem 16 If the circuit V accepts with probability at least 1 — e on some input |0£), then H 
has an eigenvalue smaller than e. If the circuit V accepts with probability less than e on all inputs 
|0£), then all eigenvalues of H are larger than 1 over a polynomial in n and L. 

The input to V consists of n\ auxiliary qubits and ni witness qubits, where n\+U2 = n. Together 
a string x of n\ bits and £ of ni bits defines an input to V. Let \(j> x Jo) be the state in template 
(n, 0) such that the input bits are set to x and Let \4> x £,t) = J^lfix^fi) for t G {1, . . . , T n }. The 
space Si e g a i is spanned by these |0)'s. 

We define 

Ke) = irp , -i Z! I ■ 

V-Ln + 1 i=0 

The following lemma established Theorem 1161 in one direction. 

Lemma 17 // there is a £ such that V accepts with probability at least 1 — e on some input |0£), 
then the smallest eigenvalue of H is at most e. 

Proof: Let \u) = |wo^). 

{v\H prop \v) = {u\Hi ega i\v) = {v\H vaUd \v) = (u\H in \u) = 0. 

For t S {0, . . . , T n — 1}, particle L is not in a G state, so the conditions of H ut are satisfied. For 
t = T n , if the system is in state \\phi ^^ t ) and the L th particle is measured, the probability that the 
outcome is Go is at most e. Therefore (^\H out \^) < e. ■ 

The final step is the following lemma: 

Lemma 18 If for all £, V accepts with probability at most e < 1/2 on input |0£), then the smallest 
eigenvalue of H in S[ ega i is at least an inverse polynomial in L. 

Proof: We will use Lemma [5j With H\ = H prop and H2 = Hi n + H out . All of Si ega i is zero on the 
other two terms, Hi ega i and pL va iid- The ground space of H prop is spanned by the \u x ^)'s, and the 
smallest non-zero eigenvalue is f2(l/(T n ) 2 ). Similarly for H2, the ground space can be expressed as 
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any linear combination of states that have G\ in particle L and do not have Q\, R\, G\ or B\ in 
particles 1 through n. Any non-zero eigenvalue is at least 1. 

Let P2 be the projection onto the ground space of H2. The cosine of the angle between any 
state \4>) and the ground state of H2 is just (0|P2 !</>)• We will show that for any |i/ Xj g), the cosine 
of the angle between \u x ^) and the ground space of H2 

Kd^Ke) < 1 - jr^i- 

Thus, for small enough e, the sine of the angle between the ground space of H\ and the ground 
space of B.2 is at least Jl(l/T n ). 

The amplitude of | 0) in \v x £) is at least l/y/T n + 1. If x 7^ 0, this state has a Qi somewhere 
in the first n particles and (4> x ,£fi\P2\4>x,£fi) = 0. Therefore, {y x ^\P2\^x,i) < 1 — y4pi ■ If x = 0, 
then we know that for all £, 

(^.TjCIGiXGiDLl^Tj^e, 
because this is the probability that the outcome on input £ is 1. This means that 
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